Nonlinear suppression of time-reversals in 'PT-symmetric optical couplers 
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We reveal a generic connection between the effect of time-reversals and nonlinear wave dynamics 
in systems with parity-time (VT) symmetry, considering a symmetric optical coupler with balanced 
gain and loss where these effects can be readily observed experimentally. We show that for intensities 
below a threshold level, the amplitudes oscillate between the waveguides, and the effects of gain and 
loss are exactly compensated after each period due to periodic time-reversals. For intensities above a 
threshold level, nonlinearity suppresses periodic time-reversals leading to the symmetry breaking and 
a sharp beam switching to the waveguide with gain. Another nontrivial consequence of linear PT- 
symmetry is that the threshold intensity remains the same when the input intensities at waveguides 
with loss and gain are exchanged. 
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I. INTRODUCTION 

The study of different physical systems exhibiting 
parity-time (VT) symmetry has attracted a lot of atten- 
tion in the past few years. The original motivation stems 
from the initial ideas of the generalization of quantum 
mechanics suggested by Bender et al. [IH^, who demon- 
strated that a wide class of Hamiltonians, even though 
non-Hermitian, can still exhibit entirely real eigenvalue 
spectra provided that they obey the so-called parity- 
time requirements or VT symmetry. A necessary (but 
not sufficient) condition for a Hamiltonian to be PT- 
symmetric is that its potential V(x) satisfies the condi- 
tion V(x) — V*(— x). It was suggested that in optics the 
refractive index modulation combined with gain and loss 
regions can play a role of complex PP-symmetric poten- 
tials [5 7]. It should be noted that the beam dynamics 
in directional couplers composed of waveguides with gain 
and loss was originally described theoretically already 
two decades ago in Ref. @, however the recently identi- 
fied analogy with VT symmetry property has stimulat- 
ing extensive theoretical and experimental HE Gil 
studies. 

As was demonstrated in the original study of direc- 
tional couplers with gain and loss [8j], such structures 
can offer benefits for all-optical switching in the non- 
linear regime, lowering the switching power and attain- 
ing sharper switching transition. Recently, these conclu- 
sions were complimented by the prediction of unidirec- 
tional switching and exact analytical solution describing 
the switching dynamics in nonlinear VT symmetric cou- 
plers [Hj]. 

In this paper, we show that the origin of nonlinear 
switching in the PT-symmetric directional couplers is re- 
lated to the effect of suppression of time-reversals. This 
conclusion is based on the symmetry analysis which is 
applicable to a broad class of nonlinear local responses, 
including in particular the cases of cubic (as considered 
in Refs. @, OH) or saturable responses. This is impor- 
tant in view of possible experimental realizations of such 
couplers in different material systems with various non- 



linear response characteristics. For example, linear PT- 
symmetric couplers have been demonstrated based on 
LiNbC>3 platform [l7|, and this material possesses pho- 
torefractive nonlinearity with saturable response. 

We show that although nonlinearity always breaks 
the PT-symmetry conditions for asymmetric wave pro- 
files even at arbitrarily small intensity levels, the ef- 
fects of gain and loss are exactly compensated and PT- 
symmetric dynamics is preserved on average due to 
periodic time-reversals, for intensities below a certain 
threshold. In contrast, for intensities above a threshold, 
nonlinear self-action suppresses time-reversals and PT- 
symmetric dynamics is broken both locally and globally, 
resulting in the asymmetric wave localization in the re- 
gion with gain. 

The paper is organized as follows. In Sec. |TT] we dis- 
cuss our model of a nonlinear PT-symmetric coupler. 
Stationary solutions of this model are found and ana- 
lyzed in Sec. IIII1 whereas Sec. IIVI provides the symmetry 
analysis underpinning general conclusions on the types 
of nonlinear dynamics on the coupler. Finally, results 
of numerical simulations which confirm and complement 
the analytical analysis are presented Sec. [V] and Sec. IVII 
concludes the paper. 



II. MODEL EQUATIONS 

We describe the propagation of waves in a PT- 
symmetric optical coupler by the equations for the mode 
amplitudes at the first and second waveguides. We use a 
set of coupled-mode equations which include additional 
terms accounting for Kerr- type nonlinearity [E [l7| : 



i— — + ipa\ + Ca 2 + G(\ai\ 2 )ai = 0, 
dz 

i ~d~z~ ~ ipa2 + + '-"(I 02 ! 2 ) 02 = 0' 



(1) 



where z is the propagation distance, a\ and a 2 are the 
mode amplitudes, p = p\ = —p 2 defines the rates of 
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FIG. 1. Scheme of nonlinear "PT-symmetric directional cou- 
pler with balanced loss in waveguide 1 and gain in waveguide 
2. The arrow indicates the propagation direction (z). 



loss in the first waveguide and gain in the second waveg- 
uide, C is the coupling coefficient between the modes of 
two waveguides, and function G characterizes the non- 
linear response. We assume with no loss of generality 
that C > 0, since for negative C it is possible to make 
the transformation a 2 — > —a 2 and C —> —C. We note 
that a recently analyzed model for Bose-Hubbard dimer 
in the mean-field approximation considered in Ref. [l8| 
may appear to be similar to the optical case, however the 
mean-field nonlinear terms are different which response 
does not depend on the total amplitude scaling in con- 
trast to Kerr-type nonlinearity in Eqs. ([1]). 

In the following, we consider the values of gain/loss co- 
efficient below the linear PT-symmetry breaking thresh- 
old [17]], p < C. In order to analyze nonlinear dynamics, 
it is convenient to represent the mode amplitudes in the 
following form, 



ai = y/l(z) cos[6>(z)] exp[+iip(z)/2] exp[i/3(z)], 



a 2 = \J I(z) sin[6>(,z)] exp[-icp(z)/2] exp[i/3(z)}, 



(2) 



where / is the total intensity, 9 and p define the relative 
intensities and phases between the two input waveguides, 
and /3 is the overall phase. After substituting Eq. @ 
into Eq. (QJ, we derive the closed system of evolution 
equations for /, 9, and p>: 



^- = -2pIcos(29), ^ = psin(26») - Csin^, 
dz dz 



dp 
dz 



(3) 



G(/cos 2 9) - G(7sin 2 



2Ccot(26>) cos ip, 



and additional equation for j3: 



df3 
dz 



1 



= - [G(I cos 2 9) /2 + G{I sin 2 9) 



C cos >p 
sin(20) ■ 



(4) 



III. STATIONARY SOLUTIONS 

Before performing the analysis of strongly nonstation- 
ary dynamics, we first analyze stationary solutions of 
Eqs. ([3]). We find that such nonlinear modes correspond 
to the phase space points (I — Iq, 9 = 7r/4, p = p>±) and 
P — P±z, where /3± = G{Iq/2) + C cos(ip±), sin(p±) = 



p/C, and cos(p±) = =f[1 - (/°/C) 2 ] 1//2 . We perform the 
linear stability analysis by considering weakly perturbed 
stationary solutions in the form: 



a 3 = exp(i/3 ± z){ v /V2exp[i(-I) J+ V±/2] 
+Uj exp(jpz) + v* exp(— ip*z)}, 



(5) 



where Uj and Vj are the perturbation amplitudes and p 
defines the growth rate. We identify non-trivial pertur- 
bation modes with eigenvalues 



p = ±2C [cos 2 ((/s) — jcos(p)] 



1/2 



(G) 



where 7 = G'(Iq/2)Io/(2C) and prime stands for the 
derivative. The fixed point in phase space is a stable cen- 
ter, lm(p) = 0, if (i) 7 < 7 cr and tp = tp_ or (ii) 7 > —% r 
and <p = tp+, where 7 cr = | cos((,c±)|. If these conditions 
are not satisfied, then the fixed point is a saddle, and a 
mode amplitude corresponding to Im(p) < grows expo- 
nentially indicating the onset of instability. 



IV. SYMMETRY PROPERTIES AND 
DYNAMICAL SCENARIOS 

In order to describe the features of nonstationary dy- 
namics, we identify an important symmetry property of 
the model equations. After performing the complex con- 
jugation of Eq. |T]) and comparing it to the original equa- 
tions, we conclude that for any solution Qj(z), 



ai{z 0+ ) = ai,(z -)e lS , a 2 {z 0+ ) = a 1 (z -)e ! 



(7) 



is also a solution of Eq. (|TJ) for arbitrary constants zq and 
5, where zq± = zq±z. This transformation represents the 
action of VT operator, where parity operator (V) corre- 
sponds to exchange of waveguide numbers and time op- 
erator (T) defines the reversal of propagation direction. 
Using notations of Eq. ([2]), we express the transformation 
in Eq. Q as 



(8) 



I ( z o+) = I(zo-), <p{z +) = p(zo-), 

9(z 0+ ) = tt/2 - e(zo-), /3(z Q+ ) = S - /3(* _). 

We notice that if zq = z m where 

9(z m ) = tt/4, (9) 

then we can choose the free parameter as S = 2/3(zq), 
and solution transforms into itself at z = z m . This 
happens because the intensity distribution is symmetric, 
|ai(zm)| 2 = |&2(z m )| 2 , anc ^ accordingly nonlinearity does 
not break the 'PT'-symmetry condition at z = z m . Since 
the original and transformed solutions satisfy the same 
evolution equation, it follows that 



I{z m +) = I(z m -), 9(z m+ ) = tt/2 - 6{z m -), 
¥>(zm+) = <p(zm-), /3(z m +) = 2f3(z m ) - P(z m -), 



(10) 
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for z rn ± = z m ± z and any z m which satisfies condition 
in Eq. ©. According to Eq. (fTU)) . the dynamics starting 
from z m in positive (+z) and negative (—z) directions 
is exactly equivalent, subject to the effective exchange of 
waveguide numbers [Eq. (J7J], and this is a nontrivial con- 
sequence of linear VI" symmetry in the nonlinear regime. 
It also follows from Eqs. ((3|) and (j9]) that 



The physical interpretation of this important result is 
that the system exhibits effective time-reversal when the 
total intensity reaches the maximum or minimum val- 
ues, where time-reversal (T) corresponds to change of the 
propagation direction (z). We use this result to reveal the 
general properties of nonstationary wave dynamics. 

First, we show that for any conditions at z = 
and considering — oo < z < +00, the dependence I{z) 
should contain a minimum (or a stationary point), where 
dl/dz = and d 2 I/dz 2 > at z min . This can be proven 
as a contradiction: (i) Depending on the sign of dl/dz 
at z = 0, we choose the positive or negative propaga- 
tion direction corresponding to decreasing intensity (if 
such direction does not exist, z = is already a mini- 
mum), (ii) As the intensity is decreasing, the effect of 
nonlincarity is diminishing (we assume the regular non- 
linear response with G(I) — > as I —> 0), and the system 
approaches linear propagation regime, (iii) However, it 
was shown that in the linear regime the intensity depen- 
dence always contains a minimum when loss coefficient 
is below the linear VT breaking threshold @, [l|| Il7j . 
Since there always exists z min , it follows from Eq. (fTUj) 
that all trajectories in phase space are symmetric with 
respect to transformation 9 — > ir/2 — 9. We note that 
for special initial conditions, at stage (ii) the solution 
may be asymptotically approaching a saddle type fixed 
point (this is not possible for the only other type of sta- 
tionary solutions in our model corresponding to a stable 
center). As we have shown above, fixed points only ap- 
pear when 9 = 7r/4, exactly fulfilling the condition of 
Eq. Therefore, such case formally corresponds to a 
limit z m i n — > +00. However, such separatrix trajectories 
approaching a saddle point are inherently unstable and 
this special case is practically unreachable and can be 
neglected. 

Second, we reveal that for arbitrary nonlinear response 
functions, all solutions belong to two classes: (i) peri- 
odic solutions, where the intensities and relative phases 
in two waveguides are exactly restored after each period 
(z — > z + z p ), or (ii) solutions where the total intensity 
grows without bound due to nonlinearly-induced symme- 
try breaking. Let us prove that solutions are periodic if 
intensity is bounded, i.e. when there exists z max where 
dl/dz = and d 2 I/dz 2 < (we again neglect the spe- 
cial case of a separatrix trajectory approaching a saddle 
point). Then, according to relation in Eq. (fTTj) . Eq. (TTU)) 
should be satisfied simultaneously for z m — z max and 



Zm — -^minj 1-6. 

I(z m in-z) = I(z min + z), I(z max -z) = I(z max +z) (12) 

Making a variable transformation z — > (z m in — z) and 
z — > (z mH - z) in the first and second relations in Eq. f| 12[) 
respectively, we obtain 

I(z) = I(2z min - z), I(z) = I(2z max - z). (13) 

Applying the second relation in Eq. (fT3| recursively after 
the first one, we find that 

I(z) =I{2z max -2z min + z). (14) 

This means that the solution is periodic, with the period 
equal to z p = 2|z max — z m - ln \. We can determine the 
location of extrema points (z max and z m - ln ) on the phase 
plane. It follows from Eq. Q and ([9]) that maxima (z max ) 
correspond to 9 = 7r/4 and ip_ < ip < (p + , and minima 
(zmin) to 9 — 7r/4 and ip < ip- or <p > tp + . 

Based on these general predictions, we can reveal a re- 
markable property. The type of nonlinear dynamics (pe- 
riodic or unbounded) remains the same if we swap the 
intensities between the two waveguides [see Eq. (fT0|) ]. In 
particular, we can couple light at the input just to the 
first waveguide with loss, or to the second waveguide with 
gain, and the type of dynamics would be the same. This 
is a counter-intuitive result, since in the first case the 
total intensity will initially decrease, whereas in the sec- 
ond case the total intensity will be growing. However, in 
both cases the type of dynamics will be determined only 
by the initial intensity level. This is a highly nontriv- 
ial consequence of linear 'PT-symmetry in the strongly 
nonlinear regime. 

V. NUMERICAL RESULTS FOR KERR 
NONLINEARITY 

We now complement the general analytical results with 
numerical examples. To be specific, we consider the Kerr- 
type nonlinear response function 

G(I) = 7 /, (15) 

where 7 > for self-focusing nonlinearity. Then, by 
introducing the transformation z — > zC and dj — > 
a j V C/7i we can scale the values of coefficients to unity, 
C = 1 and 7 = 1, and we use these values in numeri- 
cal simulations. It follows from analytical analysis pre- 
sented in Sec. Mil that the stationary point at (p+ is al- 
ways stable, whereas at if- the instability appears for 
A) > -focr = 2[1 — p 2 ] 1 / 2 . We present in Fig. [2] two exam- 
ples of system dynamics when the initial condition (at 
z = 0) is chosen in the vicinity of the unstable point, 
</?_ = 7r/6 for p = 0.5 and Iq = 2.2 > /ocr = Va. In 
both examples we choose 9 — it/ 4, such that z m = 
according to Eq. In the first example, we set the 
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FIG. 2. (Color online) System dynamics for different initial 
conditions: (a)-(d) if = tt/6 — n/20 and (e)-(h) if — ir/6 + 
7r/20. (a),(e) Trajectories in the phase plane (8, ip). Red open 
circle marks the point at z — 0, and open triangle marks the 
unstable stationary solution with ip- — tt/6. (b),(f) Intensity 
dependencies on propagation distance in the first (dotted line) 
and second (dashed) waveguides, solid line show the sum of 
individual intensities. (c),(g) and (d),(h) show the intensity 
and phase evolution along the propagation direction. For all 
the plots, p = 0.5 and I(z = 0) = 2.2. 



value of ip(z = 0) slightly larger than tp_, and therefore 
according to the general analytical results the initial con- 
dition corresponds to a maximum (z max = 0), and the 
solution should be periodic. Indeed, this is confirmed 
by numerical simulations in Fig. [^a)-(d). We see that 
the trajectory in phase space rapidly moves away from 
the unstable stationary point, but then returns to the 
initial location after a full period. Completely different 
dynamics is observed in the second example [Fig. EJe)- 
(h)], where we set the value of <p(z = 0) slightly less than 
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FIG. 3. (Color online) Regions of PT-symmetry (white shad- 
ing) and symmetry breaking with nonlinear switching (black 
shading) in the plane of initial conditions: (a)-(c) the plane 
(6, ip) for different input intensities: (a) / = O.f , (b) I = 0.8, 
and (c) I = 2.2. Stars and triangles mark stable and unsta- 
ble stationary solutions, respectively, (d) the plane (7, ip) for 
6 = 7r/4. Dashed and dotted lines mark the stationary points 
at <p±. For all the plots p = 0.5. 



ip_ . In this case, the initial condition corresponds to a 
minimum (z m i n = 0), and the solution does not have to 
be periodic. Indeed, we see that the total intensity grows 
without bound and light becomes concentrated in a sin- 
gle waveguide at \z\ — > oo. These examples illustrate two 
types of the system dynamics. 

Next, we perform comprehensive numerical studies to 
determine the dynamical regimes for arbitrary initial con- 
ditions, which may not correspond to the vicinity of sta- 
tionary points. In Fig.|3]we show the regions of the initial 
conditions corresponding to periodic (white shading) and 
unbounded (black shading) solutions, for p = 0.5. The 
plots in Figs. (3Ja-c) are shown in the plane (6, <p) for dif- 
ferent input intensities, and Fig. EJd) — in the plane (I, 
ip) for 6 — 7r/4. For the input intensity below a critical 
threshold (I < I CI ), all initial conditions correspond to 
periodic solutions preserving PT-symmetry on average. 
Note that the threshold value is lower than the instabil- 
ity threshold for the stationary mode, i.e. I cr < Ioa- As 
the input intensity exceeds the threshold, there appears 
a range of initial conditions corresponding to nonlinear 
suppression of VT symmetry reversals. On the other 
hand, according to the analytical results, the region be- 
tween the dashed and dotted lines [corresponding to the 
stationary points with tp±) in Fig. [2Id)] will always re- 
main white, since these initial conditions correspond to 
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P P 

FIG. 4. (a) Dependence of the minimal critical intensity 
required for nonlinear switching vs. the gain coefficient p. 
(b) The corresponding value of the input phase <p cr for the 
nonlinear switching. 



periodic solutions at arbitrary intensities. 

Finally, in Fig.Uwe present the dependence of the min- 
imal input intensity (J C r) on gain coefficient and the cor- 
responding phase difference ((p CI ) required for nonlinear 
switching. According to the analytical results, the mini- 
mal intensity should always correspond to the condition 



6 cr = 7r/4, and our numerical results confirm this con- 
clusion. As follows from Fig. Eta), the threshold for non- 
linear switching is drastically reduced for larger gain/loss 
coefficients, which is crucial for the practical applications 
of such couplers in all-optical signal processing. 

VI. CONCLUSIONS 

We have described analytically and demonstrated 
numerically the effect of nonlinearity-induccd VT- 
symmetry breaking in directional waveguide couplers 
with balanced gain and loss. We have revealed that time- 
reversals can support average balance between gain and 
loss despite local "PT-symmetry breaking, whereas sup- 
pression of time-reversals at stronger nonlinearities re- 
sults in switching and light concentration in a region with 
gain. Our results may offer a new insight and suggest 
different possibilities for all-optical control of the beam 
switching and amplification in nonlinear photonic struc- 
tures containing loss and gain elements 

This work has been supported by the Australian Re- 
search Council through Discovery projects. 
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